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Properties of factorial cumulant to factorial moment 



It is shown that the ratio of factorial cumulant moments to factorial moments 
for a multiplicity distribution truncated in the tail reveals oscillations in sign 
similar to those observed in experimental data. It is suggested that this effect 
be taken into account in the analysis of data in order to obtain correct physical 
information on the multiplicity distributions. 



1. Introduction 

The study of multiplicity distributions (MD's) has revealed an important empirical reg- 
ularity, as they are well described by a negative binomial distribution (NBD), in all 
reactionslH (from e+e~ annihilation to nucleus-nucleus collisions), in full phase space 
(at the highest energies, after disentangling events with different jet topology in e+e~ 
annihilation and different components in pp collisions ) and in symmetric rapidity 
intervals. On one hand, the fact that the NBD fits all available data within 10% should 
be emphasized!^; on the other hand, deviations from the fits should be carefully studied. 
One possible deviation lies in the exact shape of the tail of the distribution, which can be 
studied in general by analyzing (unnormalized) factorial moments of the distribution P„: 
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It was recently proposed to study their ratioj^: 

H, = KJF,. (3) 

which is indeed very sensitive to the properties of the tail of the distribution. The theo- 
retical basis for this analysis lies in recent results^ obtained in perturbative QCD, within 
the Modified Leading Log Approximation (MLLA), which suggest that the ratio Hq of 
eq. d^) oscillates in sign when regarded as a function of the rank q. On the contrary, 
in the case of the NBD, the ratio Hq is always positive and monotonically decreasing 
(see eq. ([T0|) below). It was pointed out in ref. |^ that many problems appear when 
one attempts to match the above mentioned QCD prediction with experimental findings; 
the most prominent ones are the inadequacy of MLLA calculations of moments of rank 
higher than the second one (higher order correlations are not under control in perturbative 
QCD), and the lack of knowledge of hadronization effects. Comparisons with experimen- 
tal data have nonetheless appeared in the recent past[^, showing qualitative agreement 
with the theoretical work. This analysis however suffers from the experimental point of 
view of one more problem that does not apply to theory 0: finite statistics does not 
allow a precise knowledge of the shape of the high multiplicity tail, because one obtains 
a truncated MD. This fact makes the straightforward experimental analysis questionable, 
as it will be shown in detail in the following section. 



2. Truncating the multiplicity distribution 

Let us consider a truncated multiplicity distribution P„, which is zero for n > uq, uq being 
the cutoff parameter. It is clear from the definition, eq. (|I]), that Fq = for q > tiq. It 
is also clear from the definition, eq. (0), that Kq is also in this case different from zero 
for all q. One can prove the following theorem about distributions and infinitely divisible 
distributions (the latter are, in our view, very important in multiparticle dynamics|§]): 

Theorem 1 : a truncated multiplicity distribution is not infinitely divisible. 

The generating function of a discrete infinitely divisible distribution can always be cast 

in the exponential form 



f{z) = e 



A[g(z)-1] 



N=0 N=0 



.m=0 



(4) 



where A is a constant, positive parameter and g{z) the generating function for the prob- 
abilities qm, so that all coefficients in the summation are non- negative. Because g{z) is a 
generating function, one can always find mo such that > 0. Then the coefficient of 
^Nmo which P/vmo IS proportioual) is always larger than or equal to e~'^A'^qJ^jj/N! for 
every > 0, therefore the distribution P„ cannot be truncated. Hence a MD which is 
truncated cannot be infinitely divisible. 

A second theorem can be proven on factorial cumulant moments' properties in a 
truncated MD: 
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Theorem 2: the factorial cumulant moments of a truncated MD are not positive definite. 
The combinants of a MD can be defined as follows 0: 



The combinants of all ranks of a MD are all non-negative if and only if the distribution is 



infinitely divisible |T0[], and since it was shown in Theorem 1 that this is not the case for 
a truncated distribution, it follows that at least one of the Cq is negative. But given the 
relation between factorial cumulant moments and combinants[^: 

oo 

^. = Ej(j-1)---(j-9 + 1)C, (6) 

3=1 

this implies that the Kg are not positive definite. Notice that since for a truncated MD 
Kg can be negative and Fg cannot, their ratio Hg can also change sign. 

In order to visually compare the effects of truncating the multiplicity distribution on 
the Hq, one can take an exact MD, truncate it and compute the resulting moments. For 
the phenomenological reasons mentioned in the introduction, we used the NBD, whose 
form and moments are well known in terms of the two parameters n and k (n is the 
average number of particles and k is linked to the dispersion D by k~^ + n'^ = D"^ /tt?): 



>NBD 



k{k + l) ...{k + n-l) f n 



n\ \n + k J \n + k 

^gik + l){k + 2)...{k + q-l) 



(7) 



F.-n' — 

' {k + l){k + 2)...{k + q-l) ^ ' 

Experimental data on charged particles MD in full phase space only present even multi- 
plicities because of charge conservation ("even-odd effect"); therefore the truncated form 
that will be discussed is 

p _ I AP^^^ if {n even) and {n < uq) ^^-^-^-^ 
1 otherwise 

Here rio is the maximum observed multiplicity and A is the normalization parameter, 
so that J2nPn = 1; the parameters n and k of the complete NBD have been chosen to 
reproduce various experimental distributions. 

In Figure |l|, Hg is plotted as a function of q for the selection of values of the parameters 
of the NBD appearing in eq. (^) listed in Table |^: the values used were taken in order 
to reproduce the MD's in e~^e~ annihilation experiments at different energies, and Uq 
was taken as the largest multiplicity in the published data for the same energy. Various 
amplitudes of oscillation are seen, depending on the parameters. Notice that the first 
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Figure 1: The ratio of factorial cumulant moments to factorial moments Hg is plotted as 
a function of the rank q for the distribution in eq. (0), for different choices of parameters, 
corresponding to a good NBD description of experimental data in e~^e~ annihilation, as 
listed in Table |T]. The lines are drawn only to guide the eye. 



Table 1: Parameters used in Figure |l]. The NBD with these values of n and k reproduces 
well the MD of e"*'e~ annihilation data in different experiments and at different cm. 
energies as listed, no is largest multiplicity in the published data. 





n 


k 


no 


Experiment 


a) 


12.9 


212.0 


28 


HRS (29 GeV) 111 


b) 


13.6 


54.0 


36 


Tasso (34 GeV) H 


c) 


15.5 


30.8 


38 


Tasso (43 GeV) m 


d) 


21.4 


24.3 


52 


Delphi (91 GeV) [0 
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Figure 2: The ratio of factorial cumulant moments to factorial moments Hg is plotted as 
a function of the rank q for the distribution in eq. (0), for different choices of parameters, 
corresponding to a good NBD description of experimental data in pp and pp collisions, as 
listed in Table |2|. The lines are drawn only to guide the eye. 



Table 2: Parameters used in Figure ^. The NBD with these values of n and k reproduces 
well the MD of pp and pp collisions data in different experiments and at different cm. 
energies as listed, no is largest multiplicity in the published data. 





n 


k 


no 


Experiment 


a) 


10.7 


11.0 


26 


ISR (30 GeV) [g 


b) 


12.2 


9.4 


32 


ISR (44 GeV) m 


c) 


13.6 


8.2 


38 


ISR (62 GeV) lH 


d) 


28.3 


3.7 


100 


UA5 (540 GeV) [g 
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Figure 3: The ratio Hq is plotted here for n = 12.2, k = 9.4 and a) no = 28, b) uq = 32, 
c) no = 36, d) no = 64. 

minimum hes around g = 5 and that the amphtude of the oscillations are comparable to 
those found in the experimental data0. 

Figure |^ is similar to Figure |I|, but the parameters are taken to reproduce MD's in pp 
and pp collider events at different energies (see Table Oscillations here are larger than 
in the previous case because a larger part of the MD is truncated away; their amplitudes 
and general shapes are still comparable to their experimental counterpart . 

The unphysical effects of the finite statistics are clearly apparent from these figures: 
this fact should be taken into consideration when analyzing the ratio Hq in experimental 
data, and these effects should be removed from the data in order to expose physical 
properties. 

In Figure ^ we show how the value of the cutoff no affects the position and amplitude 
of the oscillations: the other parameters are the same as in Figure |pD but no grows from 
28 to 64 (no = 32 in Figure]^ as in Figure ^.) In the last case no oscillations are visible 
in the figure, thus showing that such a large cutoff does not affect appreciably the ratio 
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Hq, but lower values, comparable to the experimental maximum multiplicity, give larger 
oscillations because they correspond to cutting away a larger part of the distribution. 

3. Conclusions 

We have shown that the consequences of the natural fact that the experimental multi- 
plicity distribution is truncated because of finite statistics have a large importance in the 
study of the ratio of factorial cumulant moments to factorial moments, and can mask 
physical results. An oscillatory behaviour of Hg versus q appears, which is comparable in 
size and shape to that observed in a straightforward analysis of experimental data as in 
ref. [Q. It is therefore suggested to take this effect into account in future analyses. 
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